AIAA  98-2799 

Computational  Study  of  Afterbody  Effects  for 

the  Joint  Direct  Attack  Munition  (JDAM) 

T.  L.  Donegan,  C.  T.  Bangasser,  and  J.  H.  Fox 
Sverdrup  Technology,  Inc.,  AEDC  Group 
Arnold  Engineering  Development  Center 
Arnold  Air  Force  Base,  Tennessee  37389 


19991130  089 

16th  Applied  Aerodynamics 
Conference 

June  15-18,  1998  /  Albuquerque,  NM 

For  permission  to  copy  or  republish,  contact  the  American  Institute  of  Aeronautics  and  Astronautics 
1801  Alexander  Bell  Drive,  Suite  500,  Reston,  VA  22091 


Computational  Study  of  Afterbody  Effects  for  the  Joint  Direct  Attack 

Munition  (JDAM)* 


T.  L.  Donegan,  C.  T.  Bangasser,**  and  J.  H.  Fox 
Sverdrup  Technology,  Inc.,  AEDC  Group 
Arnold  Engineering  Development  Center 
Arnold  Air  Force  Base,  TN  37389-6001 


Abstract 


Computational  analysis  of  a  JDAM  in  the  M„= 
0.9  to  1.4  flow  regime  determined  the  effect  of 
changes  to  the  store  afterbody  and  support  sting  on 
store  loads.  The  computational  results  compared 
well  to  three  sets  of  wind  tunnel  data.  Accuracy 
requirements  necessitated  both  that  the  store 
geometry  be  modeled  with  exceptional  detail,  and 
that  the  Navier-Stokes  equations  with  the  Nichols' 
k-e  turbulence  model  be  employed.  Comparisons  to 
measured  data  showed  excellent  agreement  with 
aerodynamic  load  increments  and  also  with  the 
absolute  values  of  the  loads.  Analysis  using  a  wall 
function  model  for  the  boundary  layer  also 
compared  well  with  data.  The  comparisons 
established  a  sufficient  level  of  confidence  in  the 
numerical  model  as  to  allow  its  use  in  determining 
the  forces  and  moments  free  of  support  effects  at 
the  flight  Reynolds  number. 


Nomenclature 


Aref  Reference  area 


CA 


Coefficient  of  axial  force 


CN 


Coefficient  of  normal  force 


/  Mz  y 

CLN  Coefficient  of  yawing  moment  — A — y— 

VQcxAef'ref' 

Fa  Axial  force 
Fn  Normal  force 
FY  Side  force 
/ref  Reference  length 
Mx  Rolling  moment 
My  Pitching  moment 
Mz  Yawing  moment 
M„  Free-stream  Mach  number 
pM  Free-stream  pressure 

2 

q^  Free-stream  dynamic  pressure  (l^p^V^) 
ut  Friction  velocity 

Free-stream  velocity 
yn  Normal  distance  from  wall 


CY  Coefficient  of  side  force  f — 1 

/  Mx  \ 

CLL  Coefficient  of  rolling  moment  — A — y— 

''Cl~Aref'ref-' 

(  Mv 

CiM  Coefficient  of  pitching  moment  — .  1 


y+  PwynUt/Pw 

a  Angle  of  attack,  deg 

P  Sideslip,  deg 

pw  Viscosity  at  wall 
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pw  Density  at  wall 

Free-stream  density 

Introduction 

The  harsh  environment  of  the  wind  tunnel  usu¬ 
ally  requires  relatively  large  sting  supports  for  test 
articles,  often  distorting  the  afterbody  shape.  Any 
modification  to  the  afterbody  shape  affects  the 
expected  measured  data,  necessitating  posttest 
data  analysis  to  quantify  the  effects  of  the  afterbody 
changes  and  of  the  presence  of  the  sting  itself. 

Engineers  may  determine  these  effects  experi¬ 
mentally  by  performing  an  incremental  investiga¬ 
tion  using  different  model  support  mechanisms 
such  as  struts  and  dummy  stings  to  isolate  the 
sting/afterbody  effects.  Reference  1  demonstrates 
this  technique.  This  method,  however,  assumes  lin¬ 
ear  superposition  and  requires  significant  testing 
and  model  fabrication,  escalating  program  costs. 

Rist2  established  a  computational  methodology 
that  parallels  the  experimental  approach.  He  com¬ 
puted  the  increments  for  a  sting-support  interfer¬ 
ence  problem  and  confirmed  that  the  computed 
increments  compared  well  with  measured  incre¬ 
ments.  Rist  concluded  that  one  must  accurately 
model  the  flow  physics  in  the  area  of  the  support 
geometry  and  fully  define  the  geometry  to  correctly 
determine  support  effects. 

Recently,  MDA  (McDonnell-Douglas  Aero¬ 
space)  tested  its  JDAM  (Joint  Direct  Attack  Muni¬ 
tion)  configuration  in  several  wind  tunnels  with  dif¬ 
ferent  sting  supports  and  afterbodies.  They  found 
that  different  sting/afterbody  geometries  were  nec¬ 
essary  because  of  the  different  tunnel  hardware 
arrangements  and  environments.  Because  of  the 
disparity  in  the  test  results,  MDA  requested  AEDC 
to  computationally  investigate  the  afterbody/sting 
arrangements  and  determine  the  geometric  effects 
in  the  data.  In  addition,  MDA  requested  that  scale 
effects  from  the  1  /4-scale  model  geometry  to  the 
full-scale  flight  vehicle  be  computed  and  analyzed 
as  well. 

Experimental  Data 

Three  separate  wind  tunnel  facilities  collected 
measured  data  for  the  JDAM.  The  tunnels  included 
the  McDonnell-Douglas  Poly-Sonic  Wind  Tunnel 


(PSWT),  the  Calspan/Buffalo  Transonic  Tunnel 
(TST)  and  the  Aerodynamic  Wind  T unnel  4T  of  the 
Propulsion  Wind  Tunnel  (PWT)  Facility  at  the 
Arnold  Engineering  Development  Center  (AEDC). 

The  tests  gathered  aerodynamic  load  data  at 
free-stream  Mach  numbers  of  0.4  to  1 .6  at  several 
angles  of  attack  and  sideslip.  Unit  Reynolds  num¬ 
bers  were  1 .75  million  per  foot  in  4T  and  4.0  million 
per  foot  in  the  PSWT.  Computations  used  1 .75  mil¬ 
lion  per  foot  for  the  subscale  configurations  and  20 
million  per  foot  in  the  full-scale  case. 

JDAM  Geometry 

The  JDAM  variant  for  this  analysis  is  a  modified 
Mk-84  bomb  with  a  strapped-on  guide  bomb  kit  that 
converts  it  to  a  smart  bomb.  The  kit  includes  a  sec¬ 
tioned  plate  that  conforms  to  the  Mk-84's  mid-sec¬ 
tion  surface.  Attached  to  the  plate  are  strakes 
approximately  normal  to  the  external  surface.  The 
strap-on  kit  significantly  alters  the  external  geome¬ 
try,  thus  requiring  wind  tunnel  testing  and  analysis 
to  establish  the  aerodynamic  characteristics  of  the 
modified  configuration.  Figure  1  shows  the  CFD 
surface  model  of  the  JDAM.  Modeled  with  the  Mk- 
84  bomb  are  the  strakes  of  the  strap-on  kit,  the 
attachment  lugs,  and  electrical  connector. 

The  analysis  required  four  configurations  to  be 
modeled.  These  configurations  are: 

Configuration  1:  Correct  Afterbody,  small  sting, 
1/4  scale  (CAS) 

Configuration  2:  Correct  Afterbody,  no  sting, 
full-scale,  flight  Reynolds  number  (CAN) 


Fig.  1 .  Overall  CFD  surface  definition;  correct  after¬ 
body,  small  string  case  shown. 


2 

American  Institute  of  Aeronautics  and  Astronautics 


Configuration  3:  Altered  Afterbody,  large  sting, 
1/4  scale,  span-matched  fins  (AAL2) 

Configuration  4:  Altered  Afterbody,  large  sting, 
1/4  scale,  area-matched  fins  (AAL3) 

Figures  2  to  5  depict  the  four  afterbody/sting 
configurations. 


Fig.  2.  Aft  configuration;  correct  afterbody,  small 
sting. 


Fig.  3.  Aft  configuration;  correct  afterbody,  no  sting, 
full  scale. 


Fig.  4.  Aft  configuration;  altered  afterbody,  large 
sting,  span  matched  fins. 


Fig.  5.  Aft  configuration;  altered  afterbody,  large 
sting,  area  matched  fins. 


All  four  configurations  consist  of  the  Mk-84  body 
but  have  dissimilar  afterbody/sting/fin  regions.  The 
small  sting  is  used  during  lower  dynamic  pressure 
environments,  while  larger  dynamic  pressure  con¬ 
ditions  require  the  greater  support  of  the  large 
sting.  The  large  sting  requires  the  afterbody  to  flare 
out  to  accept  the  sting,  masking  the  intersection  of 
the  fins  and  afterbody.  The  fins  may  be  modeled  by 
either  matching  the  span  and  assuming  the  fin / 
afterbody  juncture  will  be  negligible  to  the  loads  or 
by  matching  the  fin  area  and  neglecting  the  span 
and  length.  Configurations  3  and  4  represent  these 
modeling  theories. 

Computational  Model 

The  geometric  complexity  of  the  JDAM  model 
and  the  necessity  of  computing  increments  for  sev¬ 
eral  similar  configurations  required  the  use  of  the 
chimera  overset  mesh  methodology.34  The  chi¬ 
mera  methodology  allows  splitting  of  the  computa¬ 
tional  domain  into  smaller  subdomains  about  which 
body-fitted  meshes  can  be  easily  generated, 
thereby  simplifying  mesh  generation  for  complex 
geometries.  Trilinear  interpolation  of  flow-field  infor¬ 
mation  updates  the  boundary  points  of  each  subdo¬ 
main. 

The  JDAM  surface  geometry  and  near  field 
required  42  viscous  meshes,  and  the  far  field 
required  two  inviscid  meshes.  The  42  meshes  con¬ 
tained  a  total  of  2.7  million  mesh  points  for  the  1/4- 
scale  model  and  3.3  million  mesh  points  for  the  full- 
scale  model.  The  thinner  boundary  layer  of  the 
higher  Reynolds  number  flow  for  flight  required  the 
additional  points  in  the  full-scale  geometry  to  accu¬ 
rately  resolve  the  near-wall  region.  All  surface 
meshes  required  viscous  spacing  normal  to  the 
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solid  walls  with  a  dimensionless  body-normal  dis¬ 
tance  parameter  (y+)  range  between  3  and  4.  Fig¬ 
ure  6  shows  an  example  of  the  viscous  spacing,  dis¬ 
playing  the  nose  mesh  overlap  with  the  body  mesh. 

The  portions  of  the  strake  mounting  plates  that 
fit  to  the  surface  of  the  body  have  sharp,  exposed 
edges.  A  calculation  of  the  approximate  boundary- 
layer  thickness  (using  flat  plate  theory)  revealed  the 
boundary  layer  to  be  five  times  thicker  than  the 
height  of  the  plates,  indicating  the  effect  of  sharp 
edges  was  negligible.  By  filleting  the  plate  edges, 
the  plates  were  created  with  significantly  fewer 
mesh  points.  The  access  holes  and  cutouts  in  the 
mounting  plates  were  neglected  for  the  same  rea¬ 
son.  All  other  strake  mounting  hardware  (for  exam¬ 
ple,  straps,  latches,  etc.)  was  neglected  as  well. 

The  attachment  lugs  were  simplified  for  the 
computational  model  by  neglecting  the  holes  and 
modeling  the  lugs  as  solid  bodies.  Figure  7  shows 
selected  body  mesh  planes  near  the  lugs  to  dem¬ 


onstrate  the  lug  surfaces  and  the  computational 
hole  created  in  the  body  mesh.  The  electrical  con¬ 
nector  was  modeled  as  a  cylinder  that  intersects 
the  body  surface.  The  fins  were  modeled  with  a  gap 
between  the  afterbody  surface  and  the  fin  roots. 
The  fin  pivots,  however,  were  neglected. 

Wall  Function  Model 

Wall  function  analysis  of  the  JDAM  required  an 
additional  set  of  meshes  to  be  assembled  to  test  the 
recently  developed  boundary-layer  model.  The  wall 
function  model  computes  the  boundary  layer  over  a 
surface  by  relaxing  the  extremely  tight  mesh  spac¬ 
ing  requirement  normal  to  the  wall.  Using  wall  func¬ 
tions,  the  normal  wall  spacing  is  set  to  a  y+  value  of 
50  to  80,  instead  of  the  conventional  one  to  five. 
The  larger  spacing  allowed  the  flow  solution  to  con¬ 
verge  faster  because  of  larger  mesh  cells  and  fewer 
overall  mesh  points.  The  wall  function  solution  was 
approximately  two  to  four  times  faster  than  conven¬ 
tional  boundar-layer  spacing  using  the  same  solu- 
.  tion  algorithm. 


Nichols56  developed  the  wall  function 
boundary  condition  to  expedite  flow  solu¬ 
tions  for  complex  geometries.  The  boundary 
condition  models  the  inner  portion  of  the 
boundary  layer  in  conjunction  with  the  k-e 
turbulence  model  used  in  this  work.  Nichols 
has  shown  the  model  to  simulate  correctly 
the  boundary  layer  for  many  problems.  In 
this  case,  the  JDAM  was  used  as  an  exam¬ 
ple  of  a  complex  geometry  at  extreme  flow 
conditions  as  a  supplementary  test  for  the 
model,  in  addition  to  giving  credence  to  the 
remaining  solutions  in  this  work. 


Fig.  6.  Sample  of  grids;  nose  and  forebody. 


Fig.  7.  Sample  of  grids;  strake  and  lugs. 


The  wall  function  analysis  for  the  JDAM 
consisted  of  27  meshes  and  2.1  million  mesh 
points  and  was  computed  for  only  one  flow 
condition.  Compared  to  the  original  meshes, 
the  savings  in  mesh  points  and  number  of 
meshes  are  evident. 

Flow  Solution 

The  flow  solver  was  the  chimera  imple¬ 
mentation  of  a  three-dimensional  finite  differ¬ 
ence  algorithm  used  to  approximate  the 
solution  to  the  Navier-Stokes  and  Euler 
equations.  The  solver  employs  the  implicit, 
approximate  factorization  method  of  Beam 
and  Warming7.  Both  a  Baldwin-Lomax  turbu- 
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lence  model8  and  the  Nichols'  two-equation  turbu¬ 
lence  model59  are  included  in  the  flow  solver. 
Boundary  conditions  are  set  explicitly. 


dence  in  the  computations,  which  would  then  allow 
the  computational  analysis  of  full-scale  flight  Rey¬ 
nolds  number  cases. 


All  results  shown  were  computed  by  approxima¬ 
tion  of  the  Navier-Stokes  equations  with  Nichols' 
two-equation  turbulence  model.  The  wall  function 
case  employs  a  special  boundary  condition  for  the 
two-equation  turbulence  model. 

Results 

One  of  the  objectives  of  this  effort  was  the 
numerical  determination  of  the  JDAM  aerodynamic 
performance  at  flight  Reynolds  numbers.  To  estab¬ 
lish  solution  credibility  for  the  JDAM  at  flight,  where 
no  data  exist,  it  was  advantageous  to  predict  pre- 
formance  at  wind  tunnel  conditions  for  which  data 
existed.  This  required  that  numerically  determined 
aerodynamic  loads  of  three  JDAM  subscale  models 
compare  well  to  corresponding  wind  tunnel  data. 
Also,  load  differences  observed  in  the  measured 
data  due  to  afterbody  alterations  should  compare 
well  with  differences  obtained  computationally.  The 
quality  of  these  comparisons  would  provide  confi¬ 


Table  1  contains  a  list  of  CFD  solutions  com¬ 
puted  for  the  four  JDAM  configurations.  All  condi¬ 
tions  include  the  JDAM  at  a  roll  angle  of  -45  deg 
(clockwise  looking  upstream  from  lugs  up).  Wind 
tunnel  tests  previously  performed  dictated  the  con¬ 
ditions  computed. 

The  TESS10  code  enables  the  determination  of 
aerodynamic  loads  from  the  CFD  solutions.  This 
code  allows  overlapped  meshes  that  lie  on  a  sur¬ 
face  (where  pressures  may  be  double-valued)  to  be 
rectified  with  a  single  unstructured  mesh  definition, 
thus  yielding  a  unique  pressure  value  at  each  sur¬ 
face  point.  The  TESS  code  integrates  the  pres¬ 
sures  and  viscous  forces  to  obtain  the  load  compo¬ 
nents.  Figure  8  defines  the  force  and  moment  com¬ 
ponents  for  the  JDAM  geometry. 

Figures  9  to  1 1  show  comparisons  of  wind  tun¬ 
nel  data  from  the  three  tunnels  with  CFD  results  for 
the  correct  afterbody,  small  sting  configuration 


Table  1 .  Computational  Run  Matrix. 


Run 

Mach 

Config.* 

a 

p 

1 

0.9 

1 

0 

0 

2 

0.9 

1 

10 

0 

3 

1.2 

1 

0 

0 

4 

1.2 

1 

10 

0 

5 

1.4 

1 

0 

0 

6 

1.4 

1 

10 

0 

7 

0.9 

3 

10 

0 

8 

1.4 

3 

10 

0 

9 

0.9 

4 

0 

0 

10 

0.9 

4 

10 

0 

11 

1.4 

4 

0 

0 

12 

1.4 

4 

10 

0 

13 

0.9 

1 

10 

5 

14 

0.9 

4 

10 

5 

15 

0.9 

2 

0 

0 

16 

0.9 

2 

0 

5 

Configurations 

1.  Correct  Afterbody,  small  sting,  1/4 

scale  (CAS) 

2.  Correct  Afterbody,  no  sting,  full-scale, 

flight  Reynolds  number  (CAN) 

3.  Altered  Afterbody,  large  sting,  1/4 
scale,  span-matched  fins  (AAL2) 

4.  Altered  Afterbody,  large  sting,  1/4 

scale,  area-matched  fins  (AAL3) 

z 


A 


Top  View 

Fig. 8.  CFD  body  axis  coordinate  system. 
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Fig.  9.  Correct  afterbody;  small  sting,  0.9  Mach,  0  deg  |3. 
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AXIAL  FORCE 


ROLLING  MOMENT 


Angle  of  Attack  Angle  of  Attack 


Fig.  10.  Correct  afterbody;  small  sting,  1.2  Mach,  0  deg  p. 
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Angle  of  Attack  Angle  of  Attack 


Fig.  1 1 .  Correct  afterbody,  small  sting,  1 .4  Mach,  0  deg  p. 
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(CAS).  The  plots  show  good  overall  comparison  of 
CFD  solutions  with  data  at  the  three  Mach  number 
cases.  The  primary  discrepancies  also  exist  among 
the  measured  data.  These  discrepancies  occur  pri¬ 
marily  because  of  differences  in  the  Reynolds  num¬ 
bers  of  the  different  tunnel  tests.  The  computations 
adopted  the  Reynolds  number  of  the  4T  test.  In 
most  of  the  comparisons  in  Figs.  9  and  1 0,  the  com¬ 
puted  results  fall  within  the  uncertainty  of  the  4T 
instrumentation.  The  brackets  about  the  data  in  the 
figures  denote  the  uncertainty. 

In  Figs.  12  and  13,  comparisons  indicate  good 
agreement  between  CFD  and  test  data  for  the 
altered  afterbody,  large  sting,  span-matched  fin 
configuration  (AAL2).  A  comparision  of  configura¬ 
tions  shows  that  there  are  no  substantial  differ¬ 
ences  between  the  performance  of  the  AAL2  con¬ 
figuration  and  that  of  the  CAS  configuration  (com¬ 
pare  Figs.  12  and  13  with  Figs.  9  and  11).  Signifi¬ 
cant  differences  do  exist  in  pitching  moment  at  10- 
deg  angle  of  attack  between  the  AAL3  configura¬ 
tion  and  the  CAS  configuration  (compare  Figs.  14 
and  15  with  Figs.  9  and  11).  Although  there  is  a 
notable  difference  in  the  measured  pitching 
moments  between  the  configurations,  the  computa¬ 
tions  correctly  reflect  this  difference. 

By  calculating  increments  among  the  three  con¬ 
figurations,  comparisons  are  made  between  exper¬ 
imental  and  computed  increments  and  shown  in 
Figs.  16  and  17.  The  computed  differences  com¬ 
pare  very  well  to  the  test  data  differences,  espe¬ 
cially  in  the  =  1.4  cases.  However,  a  compari¬ 
sion  of  configurations  shows,  in  all  instances,  the 
large  sting,  area-matched  configuration  (AAL3), 
exhibits  the  greatest  difference  with  the  correct 
afterbody  (CAS)  configuration.  This  observation 
indicates  the  span-matched  (AAL2)  configuration 
more  closely  replicates  the  correct  afterbody.  In  this 
span-matched  case,  the  large  afterbody  flare  evi¬ 
dently  simulates  the  boundary-layer  displacement 
found  on  the  actual  afterbody  while  conserving  the 
primary  outer  fin  region.  The  area-matched  config¬ 
uration  (AAL3)  has  more  fin  extension  exposed  to 
the  flow,  producing  additional  moments  on  the 
store  body. 

Figures  18  and  19  show  comparisons  of  data 
and  computations  for  two  configurations  (CAS  and 
AAL3)  at  a  yaw  angle  of  5  degrees  and  angle  of 
attack  of  10  degrees.  The  comparisons  continue  to 


yield  excellent  agreement  of  solutions  and  data, 
even  with  a  significant  sidewash  and  total  incidence 
angle. 

Also  presented  in  Fig.  18  is  the  wall  function 
flow  solution.  The  wall  function  model  correctly  sim¬ 
ulates  the  boundary  layer,  even  at  the  relatively 
high  incidence  angle.  The  solution  produces  com¬ 
parable  results  with  the  data  and  other  presented 
solutions.  The  wall  function  model  thus  has  prom¬ 
ise  of  yielding  faster,  yet  accurate,  solutions  for 
complex  flows. 

Since  the  excellent  comparisons  between  com¬ 
putations  and  measured  data  for  all  three  configu¬ 
rations  produce  confidence  in  the  computations, 
the  free-stream  store  solution  was  computed  repre¬ 
senting  the  full-scale  geometry.  The  correct  after¬ 
body,  no-sting  case  (CAN)  configuration  represents 
the  full-scale  geometry.  Figure  20  shows  a  compar¬ 
ison  of  the  actual  afterbody  geometry  (CAN)  and 
the  correct  afterbody,  small  sting  configuration 
(CAS).  Any  differences  denote  the  effect  of  the 
sting  attachment  to  the  store.  In  this  case,  very  little 
difference  exists  except  for  the  axial  force.  The 
absence  of  the  sting  allows  closure  of  the  body  and 
the  pressure  relief  that  accompanies  it,  thus 
accounting  for  the  axial  force  difference.  This  infor¬ 
mation  indicates  the  small  sting  has  little  effect  on 
the  store  loads  at  these  conditions. 

Concluding  Remarks 

The  aerodynamic  flow  fields  about  three  sub¬ 
scale  JDAM  configurations  were  computed  to 
determine  effects  of  afterbody  changes  to  overall 
store  loads.  The  computed  results  were  compared 
to  measured  data  from  three  different  wind  tunnels. 
Overall,  the  comparisons  were  excellent,  as  both 
the  absolute  store  loads  and  differences  due  to  the 
afterbody  changes  matched  the  data  well.  Seventy 
percent  of  the  computed  forces  and  moments  fell 
within  the  stated  pre-test  wind  tunnel  uncertainties. 

The  computational  results  illustrate  the  present 
capability  at  AEDC  to  reliably  predict  the  perfor¬ 
mance  of  complex  subscale  and  full-scale  geome¬ 
tries  at  test  and  flight  Reynolds  numbers.  The  wall 
function  solution  for  the  same  configuration  addi¬ 
tionally  gives  promise  of  reliable  results  with  signif¬ 
icant  savings  of  computational  time. 
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Fig.  12.  Altered  afterbody,  large  sting,  span-matched  fins;  0.9  Mach,  0  deg  p. 
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Fig.  13.  Altered  afterbody,  large  sting,  span-matched  fins;  1.4  Mach,  0  deg  |3. 
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Fig.  15.  Altered  afterbody,  large  sting,  area-matched  fin;  1.4  Mach,  0  deg  p. 
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Fig.  17.  Calculated  increments;  1.4  Mach,  0  deg  p. 
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Fig.  18.  Correct  afterbody,  small  sting;  0.9  Mach,  5  deg  (3. 
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Fig.  19.  Altered  afterbody,  large  sting,  area-matched  fins;  0.9  Mach,  5  deg  p. 


17 

American  Institute  of  Aeronautics  and  Astronautics 


References 


1.  Whitby,  D.  G.  "Wind  Tunnel  Support  System 
Effects  on  a  Fighter  Aircraft  Model  at  Mach  Num¬ 
bers  from  0.6  to  2.0."  AEDC-TR-89-4  (AD- 
A210614),  July  1989. 

2.  Rist,  M.  J.  "Computational  Investigation  of 
Wind  Tunnel  Support  Interference."  AIAA  94-2602. 

3.  Benek,  J.  A.,  Steger,  J.  L.,  Dougherty,  F.  C., 
and  Buning,  P.  G.  "Chimera:  A  Grid-Embedding 
Technique."  AEDC-TR-85-64  (AD-A1 67466), 
December  1985. 

4.  Benek,  J.  A.,  Donegan,  T.  L.,  and  Suhs,  N.  E. 
"Extended  Chimera  Grid  Embedding  Scheme  with 
Application  to  Viscous  Flows."  AIAA  87-1126CP, 
June  1987. 

5.  Nichols,  R.  H.,  "Development  and  Validation 
of  a  Two-Equation  Turbulence  Model  with  Wall 
Functions  for  Compressible  Flow."  AIAA  96-2385, 
June  1996. 


6.  Suhs,  N.  E.,  Nichols,  R.  H.,  and  Denny,  A.  G., 
"Unsteady  Viscous  Flow  Computations  Using  a 
Two-Equation  Turbulence  Model  with  Wall  Func¬ 
tions."  AIAA  96-2430,  July  1996. 

7.  Beam,  R.,  and  Warming,  R.  F.,  "An  Implicit 
Finite-Difference  Algorithm  for  Hyperbolic  Systems 
in  Conservation-Law  Form."  Journal  of  Computa¬ 
tional  Physics,  Vol.  22,  No.1,  September  1976,  pp. 
87-110. 

8.  Baldwin,  B.  S.  and  Lomax,  H.  “Thin  Layer 
Approximation  and  Algebraic  Model  for  Separated 
Turbulent  Flows."  AIAA  78-257,  January  1978. 

9.  Nichols,  R.  H.  "A  Two-Equation  Model  for 
Compressible  Flows."  AIAA  90-0494,  January 
1990. 

10.  Dietz,  W.  E.  "General  Approach  to  Calculat¬ 
ing  Forces  and  Moments  on  Overset  Grid  Configu¬ 
rations."  Proceedings  of  the  Second  Overset  Com¬ 
posite  Grid  and  Solution  Technology  Symposium, 
October  1 994. 


19 

American  Institute  of  Aeronautics  and  Astronautics 


